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Nonlinear g-voter model with inflexible zealots 


Mauro Mobilia^ 

^Department of Applied Mathematics, School of Mathematics, University of Leeds, Leeds LS2 9JT, 

We study the dynamics of the nonlinear g-voter model with inflexible zealots in a finite well-mixed 
population. In this system, each individual supports one of two parties and is either a susceptible 
voter or an inflexible zealot. At each time step, a susceptible adopts the opinion of a neighbor if 
this belongs to a group of g > 2 neighbors all in the same state, whereas inflexible zealots never 
change their opinion. In the presence of zealots of both parties the model is characterized by a 
fluctuating stationary state and, below a zealotry density threshold, the distribution of opinions is 
bimodal. After a characteristic time, most susceptibles become supporters of the party having more 
zealots and the opinion distribution is asymmetric. When the number of zealots of both parties is 
the same, the opinion distribution is symmetric and, in the long run, susceptibles endlessly swing 
from the state where they all support one party to the opposite state. Above the zealotry density 
threshold, when there is an unequal number of zealots of each type, the probability distribution is 
single-peaked and non-Gaussian. These properties are investigated analytically and with stochastic 
simulations. We also study the mean time to reach a consensus when zealots support only one party. 


PACS numbers: 89.75.-k, 02.50.-r, 05.40.-a, 89.65.-s 

I. INTRODUCTION 


The voter model (VM) [^j is one of the simplest and 
most influential examples of individual-based systems ex¬ 
hibiting collective behavior. The VM has been used as 
a paradigm for the dynamics of opinion in socially inter¬ 
acting populations, see e.g. iJ and references therein. 
The classical, or linear, VM is closely related to the Ising 
model Q and describes how consensus results from the 
interactions between neighboring agents endowed with a 
discrete set of states (“opinions”). While the VM is one 
of the rare exactly solvable models in non-equilibrium 
statistical physics, it relies on oversimplified assumptions 
such as perfect conformity and lack of self-confidence of 
all voters. This is clearly unrealistic as it is recognized 
that members of a society respond differently to stimuli: 
Many exhibit conformity while some show independence, 
and this influences the underlying social dynamics iil- 
In order to mimic the dynamics of socially interacting 
agents with different levels of confidence, this author in¬ 
troduced “zealots” in the VM [p|-^. Originally zealots 
were agents favoring one opinion ^ . The case of in¬ 
flexible zealots whose state never changes was then also 
studied and the influence of committed and/or inde¬ 
pendent individuals was considered in various models of 
opinion and social dynamics OEl- Recently, authors 
have investigated the effect of zealots in naming and co¬ 
operation games, and even in theoretical ecology Oil- 
In recent years, many versions of the VM have been 
proposed [^. A particularly interesting variant of the 
VM is the two-state nonlinear g-voter model (qVM) intro¬ 
duced in iQ. In this model q randomly picked neighbors 
may influence a voter to change its opinion. When q = 2, 


the gVM is closely related to the Sznajd model |I5-I7 
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and to that of Ref. . The properties of the qYM have 
received much attention and there is a debate on the ex¬ 
pression of the exit probability in one dimension [p^[-pO| . 

Here, we investigate a generalization of the nonlin¬ 
ear gVM, with g > 2, in which a well-mixed popula¬ 
tion consists of inflexible zealots and susceptible voters 
influenced by their neighbors. As a motivation, this par¬ 
simonious model allows to capture three important con¬ 
cepts of social psychology [|| and sociology [||: (i) con¬ 
formity/imitation is an important social mechanism for 
collective actions; (ii) group pressure is known to influ¬ 
ence the degree of conformity, especially when a group 
size threshold is reached |Q ; (iii) the degree of conformity 
can be radically altered by the presence of some individ¬ 
uals that are capable of resisting group pressure i I- 
Here, the gVM mimics the process of conformity by imi¬ 
tation with group-size threshold, whereas zealots are in¬ 
dependent agents that resist social pressure and can thus 
prevent to reach unanimity. 

In this work, we study the fluctuation-driven dynam¬ 
ics of the two-state gVM with zealots in finite well-mixed 
populations and shed light on the deviations from the 
mean field description and from the linear case (q = I). 
We find that below a zealotry density threshold the prob¬ 
ability distribution is bimodal instead of Gaussian and, 
after a characteristic time, most susceptibles become sup¬ 
porters of the party having more zealots. When both 
parties have the same small number of zealots, suscepti¬ 
bles endlessly swing from the state where they all support 
one party to the other with a mean switching time that 
approximately grows exponentially with the population 
size. 

In the next section we introduce the model. Sections 
HI and IV are dedicated to the mean field description 
and to the model’s stationary probability distribution. 
In Secs. V and VI we discuss the long-time dynamics 
and the mean consensus time when there is one type 
of zealots. We summarize our findings and conclude in 
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Sec. VII. 


II. THE g-VOTER MODEL WITH ZEALOTS 

We consider a population of N voters that can support 
one of two parties, either A or B, and therefore be in two 
states. Supporters of party A are in state +1, and those 
supporting party B are in state —1. Among the voters, 
a fixed number of them are “inflexible zealots” while the 
others are “susceptibles”. Here, zealots are individuals 
that never change opinion: they permanently support 
either party A (A-zealots) or party B (B-zealots). Sus¬ 
ceptible voters can change their opinion under the pres¬ 
sure of a group of neighbors. The population thus con¬ 
sists of a number of A-zealots (pinned in state -1-1) 
and Z_ of B-zealots (pinned in state —1), and a total of 
S = N — — Z- susceptibles agents, of which n are A- 

susceptibles (non-zealot voters in state -1-1) and S — n are 
B-susceptibles (non-zealot voters in state —1). The frac¬ 
tion, or density, of susceptibles in the entire population 
remains constant and is given by s = S/N. For simplic¬ 
ity we assume that all agents have the same persuasion 
strength. 

At each time step, a susceptible voter consults a group 
of q neighbors (with q > 1) and, if there is consensus in 
the group, the voter is persuaded to adopt the group’s 
state with rate 1 0- The dynamics is a generaliza¬ 
tion of the nonlinear gVM with a finite density of 
zealots 0, and consists of the following steps: 

1. Pick a random voter. If this voter is a zealot noth¬ 
ing happens. 

2. If the picked voter is a susceptible, then pick a 
group of q neighbors (for the sake of simplicity rep¬ 
etition is allowed, as in Refs. [0, 0). If all g neigh¬ 
bors are in the same state, the selected voter also 
adopts that state. Nothing happens in the update 
if there is no consensus among the q neighbors [0 , 
or if the voter and its q neighbors are already in 
the same state. 

3. Repeat the above steps ad infinitum or until con¬ 
sensus is reached. 

The case <7 = 1 corresponds to the classical (linear) voter 
model &I-0 , and we therefore focus on g > 2. 

For the sake of simplicity, we investigate this model on 
a complete graph (well-mixed population of size N). The 
state of the population is characterized by the the prob¬ 
ability P{n,t) that the number of A-susceptibles at time 
t is n. This probability obeys the master equation 0| 

= r+_iP„_i(t) + T-+iP„+i(t) 

- (r+ + T-)p„(t). (1) 

The first line accounts for processes in which the number 
of A-susceptibles after the event equals n, while the sec¬ 
ond term accounts for the complementary loss processes 
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FIG. 1: Number n{t) of A-susceptibles vs. time in two 
sample realizations (black and gray) at low zealotry, {Z+ -|- 
Z-)I2N < Zc, with initial condition n(0) = S/2, see text. 
Here, g = 2, A = 200 and S = 120. (a) Symmetric zealotry 
with Z+ = Z- = 40, Zc = 0.25: n{t) continuously fluctu¬ 
ates and suddenly switches from n « S to n « 0 and vice- 
versa. (b) Asymmetric zealotry with Z+ — 43 and Z- = 37, 
Zc = 0.2022: After a transient, n{t) fluctuates around a value 
corresponding to a majority of A-susceptibles, see text. 


where n —> n ± 1. Here, represent the rates at which 
transitions occur and are given by 


/S — n\ /n + Z+ 
V N ) \ N-1 


\ _ n f S + Z_-n 

’ ” N \ N-1 


( 2 ) 


When there are zealots of both types = T~^q = 0 
and the system has reflective boundaries at n = 0 and 
n = S. When there are only A-zealots, Z- = 0 and 
Z+ = N( > 0, with C being the density of A-zealots, 
then n = S = N{1 — (/) is an absorbing boundary with 
= R while n = 0 is reflective. The birth-and-death 
process (|^ is here simulated with the Gillespie algorithm 
0 upon rescaling time in Eq. (|^) as t —>■ t/N. 

A quantity of particular interest is the magnetization 
TO = [n -I- Z+ — {S — n) — ZN\/N = [2n — S — 
Z-]/N that gives the population’s average opinion or, 
equivalently here, the opinion of a random voter [0. We 
have TO = TOinax = {S+Zj^ — Z-)/N when all susceptibles 
are state -1-1 (all A-susceptibles) and to = TOmin = (—5'-!- 
Z+ — Z-)/N when all susceptibles are in state —1 (all 
B-susceptibles), with TOmin < ™ < ITlmax- 

To gain an intuitive understanding of the gVM dynam¬ 
ics, it is useful to consider the evolution of n(t) in typical 
sample realizations, as those in Fig. where we illustrate 
the dynamics at low zealotry. In Fig. we notice two 
distinct regimes and different time-scales. In the case of 
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FIG. 2: Schematic of the mean field dynamics when z+ = 
Z- = z and x* = s/2 (a,b); when z± = {1 ± S)z > 0 with 
0 < 5 < 1 (c,d); and when 2 + = ^ > 0,= 0 (e,f). (□) 
and (•) indicate stable and unstable fixed points, respectively. 
Panels in the left column correspond to low zealotry z < Zc 
(a,c) and C < Cc (e); those in the right column correspond to 
high zealotry z > Zc (b,d) and C > Cc (f): see text. 


symmetric low zealotry (Z+ = Z-), the number of sus- 
ceptibles first approaches either the state n ~ 0 (all B- 
susceptibles) or n « S' (all A-susceptibles). After a char¬ 
acteristic time (see Sec. V.A), all susceptibles suddenly 
start switching from one state to the other, see Fig. (a). 
A similar feature has been observed in the Sznajd model 
(q = 2) with anticonformity 0 - When > Z- > 0, 
the majority of susceptibles become A-supporters after a 
typical time (see Sec. V.B). The fluctuations in the num¬ 
ber of A-susceptibles then grow endlessly, see Fig. 0 (b). 
An important aspect of this work is to analyze how de¬ 
mographic fluctuations arising in finite populations alter 
the mean held predictions. In Section V the phenomena 
illustrated by Fig. 0 are studied in large-but-hnite popu¬ 
lations, and we show that these phenomena are beyond 
the reach of the next section’s mean held analysis. 


III. MEAN FIELD DESCRIPTION 

For further reference, it is useful to consider the mean 
held (MF) limit of an inhnitely large population, A —>■ oo. 
In such a setting, demographic huctuations are negligible 
and the rates ( 0 ) can be written in terms of the density 
X = n/N of A-susceptibles, and the densities z± = Z±/N 
of zealots of each type: T+ —)> T'^{x) = (s — x)(x + 2 +)'^ 
and T~ —>■ T~(x) = x(s + Z- — x)'^. The MF dynamics 
is described by the rate equation obtained by averaging 
n/N from Eq. (0) (and rescaling time a.s Nt ^ t) 

X = T+{x)-T-{x) 

= {s — x){x + z+y — x{s — X + Z-Y, (3) 

where the dot denotes the time derivative and s = S/N. 

In the absence of zealotry {z± = 0, s = 1), Eq. (0) 
has two stable absorbing hxed points, x = 0 (all B- 
supporters) and x = 1 (all A-supporters) correspond¬ 
ing to consensus with either A or B party, separated by 
an unstable hxed point x = 1/2 (mixture of A- and B- 
voters) 1^. It is worth noting that the dynamics of the 


gVM without zealots ceases when a consensus is reached 
and this happens in a hnite time when the population 
size is hnite UMl- However, in the presence of zealots 
supporting both parties, the population composition end¬ 
lessly huctuates 00, see, e.g., Fig. 0. 

In the presence of zealotry, the interior hxed points of 
Eq. (0) satisfy T+(x) = T“(x), which leads to 


/s — x\ / X-I-Z+ 

\ X ) \s — X + Z- 


9 

= I. 


(4) 


Depending on the values of z± and q, this equation has 
either three physical roots, or a single physical solution. 


A. The symmetric case z+ = Z- = z 


When the density of zealots of both types is identical, 
z+ = Z- = z and s = I — 2^ with 0 < z < 1/2, Eq. 0) 
becomes 


X = (I — 2z — x)(x -I- zy — x(I — z — x)'^. 


that is characterized by a hxed point x* = s/2. When z 
is sufficiently low, Eq. (0) has two further hxed points: 
x/^_ and x*_ = s — x!j_. The analysis for arbitrary 9 > I is 
unwieldy, but insight can be gained by focusing on q = 2 
and 9 = 3, for which 

* _ / ± y i - 4z ) (<7 = 2) 

~ i ^ (s ± (9 = 3). 

We readily verify that x^ are both stable when z < Zc{q), 
with Zc(2) = 1/4 and Zc(3) = 1/3. When z > Zc{q), the 
hxed points Xj- are unphysical and x* = s /2 is stable. 
This picture holds for arbitrary hnite value oi q > 1: 
Xj- are stable and the MF dynamics is characterized by 
bistability below a critical zealotry density Zc{q), while 
X* = s/2 is unstable when z < Zc and stable when z > Zc, 
see Fig. 0(a,b). By determining when Eq. 0) has three 
physical roots, we have found the critical zealotry density 


r 1/4 (q = 2) 

I 1/3 (q = 3) 

1 3/8 {q = 4) 

[ 2/5 {q = 5), 


(5) 


while Zc(l) = 0 since in the linear VM Eq. (0) has always 
one single stable hxed point ||^. Hence, the value of 
Zc increases with q, while the values of x\ and x*_ get 
closer to the values 0 (all B-susceptibles) and s (all A- 
susceptibles) as q increases with z kept hxed. 

In this MF picture, the population’s average opinion 
given by the magnetization m{t) = 2 x(t) — s undergoes 
a supercritical pitchfork bifurcation at z = Zc |^: At 
t —>■ 00 , the critical value Zc separates an ordered phase 
(z < Zc), where a majority of susceptibles supports one 
party, from a disordered phase (z > Zc) in which each 
party is supported by half of the susceptibles, see Fig. 0 
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depends on the initial condition and is m(oo) = = 

2(a:+ + Sz) — s if to(0) > 2{x* + Sz) — s and m{oo) = 
m*_ = 2{x*_ + 5z) — s if m(0) < 2{x* + Sz) — s. When 
z > Zc the stationary MF magnetization is m{oo) = m\. 


C. The absorbing case z+ = C,, Z- = Q 

When there are only A-zealots, = C > 0 and Z- = 0, 
Eq. (^) becomes 

i = (1 - c - a;) [(a; + C)'^ - a:(l - c - , 


FIG. 3: {Color Online) Critical value of Zc as a function S = 
{z+ — z-)l2z for g = 2 (solid), q — 3 (dashed), and g = 4 
(dash-dotted) in the case of asymmetric zealotry. There is 
bistability where « = {z+ — Z-)/{25) < Zc, see text. 


(a,b). The stationary MF magnetization thus depends on 
the initial condition: when 2 < Zc, m(oo) = m* = 2x’^ — s 
if m(0) > 0 and m(oo) = —m* if m(0) < 0, while the 
magnetization vanishes when z > Zc (or if m(0) = 0). 
Using Eqs. (||) and (||), it can be directly checked that 
just below the critical zealotry density, i.e. for z < Zc, the 
stationary magnetization is characterized by the scaling 
relationship m(oo) oc m* ^ y/zc — z. 


B. The asymmetric case z± = (1 ± 5)2 

When the number of A-zealots exceeds that of B- 
zealots, with z+ > z_ > 0, it is convenient to use the 
parametrization 


z± = {l±S)z, (6) 

where 6 = (z_|_ — z_)/2z quantifies the zealotry asym¬ 
metry. With Eq. (^), we still have s = 1 — 2z with 
0 < z < 1/2 and Eq. (^) becomes 

X = (1 — 2z — a:)[a: -I- (1 + S)z]‘^ — x[l — (1 + S)z — x]‘^. 

This rate equation is also characterized by bistability at 
low zealotry, with two stable fixed points a;5i separated by 
an unstable fixed point x*, and by the sole stable fixed 
point x\ at higher zealotry, see Fig. ^c,d). By deter¬ 
mining when Eq. (^ has three physical fixed points, we 
have determined the critical density of zealotry Zc(g,5), 
see Fig. |3[ At fixed q and S, the fixed points x\ are sta¬ 
ble when z < Zc while only a;)j_ is stable when z > Zc. 
We have found that Zc decreases with S (at fixed q) and 
increases with q (at fixed <5). 

In this MF picture, the opinion of a random individual 
is given by the magnetization m{t) = 2{x + Sz) — s. The 
critical zealotry density Zc separates a bistable phase (z < 
Zc) from a phase where most susceptibles support the 
party having more zealots, see Fig. ^ (c,d). Hence, the 
stationary MF magnetization at low zealotry (z < Zc) 


and has an absorbing fixed point x* = 1 — (/. Below a 
critical zealotry density Cc{<l), this rate equation admits 
two other fixed points: x^, that is stable, and x* that is 
unstable and separates x* and x^, see Fig. ||(e,f). When 
C > Cc(q), the absorbing state cc* = 1 — ^ is the only fixed 
point. For q = 2 and g = 3, we explicitly find 


i (l - 3C - ^/l-(6-C)C) 

l-2C(l-^C)-vT^ 

2f2+C) 


(q = 2) 
(q = 3) 


(7) 


and 


i(l-3C+v'l-(6-C)C) 

2(2-HC) 


(q = 2) 
(q = 3) 


( 8 ) 


From these expressions, and more generally by determin¬ 
ing when Eq. (|^) has three physical fixed points, we have 
found the critical zealotry density in the absorbing case: 


Cc(q) 


3 - 2 V 2 (q = 2) 

1/4 (q = 3) 

0.295 (q = 4) 

0.326 (q = 5) 


(9) 


We thus distinguish two regimes: 

(i) When ( < (c(q) both a:*j, are stable and the dy¬ 
namics crucially depends on the initial density Xg of A- 
susceptibles: If xq > x*, the final state is the consensus 
with party A; whereas the steady state consists of a vast 
majority of B-party voters when xg < x*. In Sec. VI, 
we show that random fluctuations drastically alter this 
picture: In a finite population, xl is a metastable state 
when ( < Cc{q) and xg < x*, and we shall see that the 
A-consensus is reached after a very long transient that 
scales exponentially with the population size. 

(ii) When ( > Cc(g), as well as when C = Cc and xg > 
X*, the absorbing state is rapidly reached. 


IV. STATIONARY PROBABILITY 
DISTRIBUTION 

In this section, we compute the stationary probability 
distribution (SPD) of the gVM with zealotry when there 
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is no absorbing state, and show that it shape generally 
differs from the Gaussian-like distribution obtained in the 
linear VM with zealots [p^ . 

The SPD P* = lim(_>.oo Pn (t) obeys the following sta¬ 
tionary master equation, obtained from Eq. (|^): 

t+_,p:_,+t-^,p:+, - ( r ++ t-)p: = o . 

The exact SPD is uniquely obtained by iterating the de¬ 
tailed balance relation T^_iP*_i = T~ P* |^, yielding 

n —1 

Pn = PS 



— p* 

— -in 


i-0 
n —1 

n 

3=0 


S-j \ ( J + Z+ 

3 + \) \s + Z_-j-l 


( 10 ) 


where the normalization X)n=o Pn gives Pq = 1/[1 -|- 
ELi /T-+^)] and P* = 1 - P* - yIzI PS- 

Since n = N[{m -I- s)/2 — 5z], the stationary magneti¬ 
zation distribution QSi has the same shape as PS, with 


Qm PN[(m+s)/2—5z\ 

N[{m-\-s) l2 — 5z] — l 


_ p* 

— -in 


n 

1=0 


j + ^4 


S-j 

j + lj \S + Z.-j-1 


In large populations, a useful approximation of ( 0 ) 
is obtained by writing P* = Pq* exp 'I'l) 

'i’j = \n{T^/T~^-^), and by using Euler-MacLaurin for¬ 
mula E”=o ^1 = Io~^ '^3 + (^0 + 5'n-i)/2, where we 

have neglected higher order terms 

When N ^ 1, it is useful to work in the continuum 
limit with the rates —>■ T^{x), as in Sec. III. By 

introducing 


vl/(cr) =ln[T+(cr)/T-(a;)], 


( 12 ) 


CL 

CO 


( 11 ) 



FIG. 4: Rescaled SPD at low zealotry with = Z- = z < Zc 
in semi-log scale, (a) sPS vs. n/s from Eq. ( |lCl| ) for different 
values of q and z, with N = 200. Here, q = 2, z = 0.2 (O); 
q = 3, z = 0.22 (A); q = 3, z = 0.2 (o) and g = 4, ^ = 0.2 (□). 
Lower inset: Similar; sPS vs. n/s from stochastic simulations. 
Upper inset: J// dy vs. x/s for z = 0.2 and, from top 

to bottom, g = 1 (dashed), g = 2 (black), g = 3 (gray) and 
g = 4 (light gray), (b) SP/ vs. n/S from Eq. ( |lo[ ) for g = 2 
(diamonds) and g = 3 (circles), and for different values of the 
population size N at low zealotry. Here, 2 = 0.2 and N = 200 
(open symbols), N = 300 (gray-filled symbols) and N — 400 
(black-filled symbols). Not shown in panels (a) and (b) is the 
range where P/ < 10~® (where P/ < 10~® in the lower inset). 


we have fg dx to leading order in N. 

Hence, the leading contribution to the SPD when iV >> I 
is 


P* 


PS exp (^N ^ dc(j/) dy^ = P*(x). (13) 


The local extrema of P*{x) satisfy 'F(a;) =0, see (1^), 
and thus coincide with the fixed points of Eq. (||). As 
a consequence, in large populations P* is either charac¬ 
terized by a single peak at n* = Nx* when z > Zc, or 
has two peaks at the metastable states n’S- = Nx'^ when 
z < Zc- In this case, there is bistability and the ampli¬ 
tudes of the peaks at are in the ratio {N :s> 1) 




'l'(y) dy 


(14) 


The integrals in Eqs. (^ and ( p^ can be computed, but 
their expressions are unenlightening. Here, we infer the 
properties of P* ~ P*ix) and from those of 4'(a:). 


A. Stationary probability distribution in the 
symmetric case 


In the symmetric case, = z_ = z, Eq. becomes 


fE'(a;) = In 


I — (a; -I- 2z) 


X + z 


I — {x + z) 


and has the symmetry 'I'(x) = —'I'(s —a:). We distinguish 
the cases of low and high zealotry density: 

(i) When z < Zc{q), the fixed points x* and x*^ of 
Eq. (^ are also the roots of 4'(a;). Hence, when N ^ 
I, P* = Ps-n ^ P*i^) oc is a symmetric 

bimodal SPD characterized by two peaks at n = nj_. As 
a consequence, QSi = Q-m is an even function. 

In Figure ^ (a), we show the exact SPD for g = 2 — 4 
characterized by two peaks of same intensity at n = n^. 
and a local minimum at n* = 3/2. We remark that when 






















6 



FIG. 5: sPn vs. n/s for g = 2 — 4 at high zealotry 2 : > Zc- 
Here, N = 200 and {q,z) = (2,0.4) (O), (3,0.4) (o), (3,0.375) 
(A), (4, 0.4) (□). The SPDs have a single peak a,t njs = N/2 
and width broadens when q and 1/z are increased at fixed N. 
Inset: SP^ vs. n/S for q = 2,z = 0.4 and different values of 
N. Here, N = 200 (O) and N = 600 (♦). 


q is increased, the SPD vanishes dramatically away from 
the peaks. In fact, since '^/{y)dy is close to zero or 
negative on x*_ <C x <C see Fig. ^ (a, upper inset), 
P*, vanishes exponentially with N and when q 

is increased. Fig. ||(a) shows that the SPD steepens and 
its peaks are more pronounced when q and 1 /z are in¬ 
creased and N is kept fixed. We have also obtained the 
(quasi-)SPD from stochastic simulations, see Fig. ^ (a, 
lower inset), by averaging over 25,000 realizations af¬ 
ter 40,000 simulation steps. While unavoidably more 
noisy, the simulation results reproduce the predictions of 
Eq. (|l^). Fig. ^(b) shows how SP* scales with n/S = x/s 
for different population sizes, and we notice that the main 
influence of raising N is to concentrate the probability 
density SP* around the peaks whose location are essen¬ 
tially unaffected by N (when iV ^ 1). In Fig. H, we also 
notice that the symmetric peaks are clearly identifiable 
when q = 2 and q = 3, but almost coincide with n = 0 
and n = S ioT q = 4. This is because Xj- approach the 
values X = 0, s when q is increased. 


(ii) When 2 ; > Zc, the only physical root of 41(x) is 
X* = s/2, as in the classical voter model [|^. Hence, 
P* = p*_^ ^ p*{x) oc has a single maxi¬ 

mum at X = s/2 when N ^ 1. The resulting symmetric 
Gaussian-like distribution centered at n* = S/2 when 
N ^ 1, see Fig. ||, is very similar to the SPD obtained in 
the classical voter model with zealots Fig. ||(inset) 
illustrates that the probability density steepens around 
s/2 when the population size is increased. 



n/s 



FIG. 6: Rescaled SPD under asymmetric zealotry Z± = 
N{1 ± 5)z. (a) sP/ vs. n/s from Eq. (^^ at low zealotry 
{z < Zc). Here, N = 200 and {q, Z+, Z-) = (2,41,39) (O), 
(3,51,49) (o), (3,53,47) (A), (4,51,49) (□).For g = 4, the 
range where P/ < 10“^^ is not shown. Inset: SP/ vs. n/S for 
g = 2 and {N, Z+, Z-) = (200, 41, 39) (open O), (400, 82, 78) 
(gray-filled symbols), (b) Left-skewed rescaled SPD at high 
zealotry, with a single peak at nlj_. Here, N = 200 and 
{q,Z+,Z.) = (2,67,63) (O), (3,72,68) (o), (3,74,66) (A), 
(4,82,78) (□), see text. 


B. Stationary probability distribution in the 
asymmetric case 


In the asymmetric case, with zealot densities z± = 
(1 ± < 5 ) 2 : and 5 > 0, Eq. ( p^ is 


4'(x) = In 


^1 — (x -I- 2z)^ 


X -b z(l -b 6 ) 

1 — z{l -b i5) — X 


91 


and has either three or one physical roots: 

(i) At fixed q and 5, when z < Zc(g, <5), the fixed points 
X* and Xj- of Eq. (||) are the physical roots of 'I'(x). 
Since P* ^ P*{x) oc e^-^o '^iv)d-v when N ^ 1, the 
SPD is again a bimodal distribution peaked at nj_. How¬ 
ever, x*c^ has a greater basin of attraction than x*_ and 

'F(j/) dy >0. As a consequence, the SPD is asym¬ 
metric, with the peak at n!/ being much stronger than 
the one at The ratio of the peaks is given by (0): 
which shows that the asymmetry of P* grows exponen¬ 
tially with N and increases with q, see Fig. ^(a). While 
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an asymmetry in the zealotry in the linear VM does not 
significantly affect the form of the SPD [Q , we here find 
that in the gVM even a small bias in the zealotry drasti¬ 
cally changes the shape of the SPD and leads to marked 
dominance of the party with more zealots. 

In Fig. ^ (a), we report the exact SPD for g = 2 —4 and 
illustrate its asymmetric bimodal nature, with marked 
peaks of different intensities at n’^. We notice that 
the asymmetry in the peaks intensity, given by & , is 
stronger when we increase q and — Z- (x 6z. As in the 
symmetric case, the SPD decays dramatically away from 
the peaks and P*. vanishes with N 1 and when 

q is increased. In Fig. || (a, inset) we show that the SPD 
remains bimodal when the population size is increased, 
and the main influence of raising N is to concentrate the 
probability density near its peak at = n*^_/N (when 
iV» 1). 

(ii) At fixed q and S, when z > Zc, the only real root 
of 'I'(a;) is a; = This lies closer to a; = s than to 
a; = 0, and hence f^^(y)dy is an asymmetric function 
with a single maximum at a; = Therefore, in large 
populations P* is an asymmetric left-skewed SPD with 
a single peak at n = n^, as shown in Fig. ib) where we 
see that the SPD broadens when q is increased and that 
it steepens when — Z- oc Sz is increased. As above, 
the probability density steepens around when N is 
increased. 


V. FLUCTUATION-DRIVEN DYNAMICS AT 
LOW ZEALOTRY 

We now study how a small non-zero density of zealots 
of both parties (0 < z < Zc) affects the gVM long-time 
dynamics. We show that, after a typically long tran¬ 
sient, all susceptibles voters switch allegiance from the 
state n = 0 (all B-susceptibles) to state n = S (all A- 
susceptibles) in a typical switching time. In the symmet¬ 
ric case, there is “swing-state dynamics” with all suscep¬ 
tibles endlessly swinging allegiance. In the asymmetric 
case where party A has more zealots than party B, the 
dynamics is characterized by various time-scales and by 
growing fluctuations around the metastable state . Be¬ 
low, we show that the long-time gVM dynamics is driven 
by fluctuations and characterized by a mean switching 
time that scales (approximately) exponentially with the 
system size N in large-but-finite populations. 


A. Swing-state dynamics and switching time in the 
case of symmetric zealotry 

As illustrated in Fig. |^(a), the long-time dynamics 
in the symmetric case is characterized by the contin¬ 
uous swinging from states n ~ 0 to n si S' and vice 
versa. When Z+ = Z-, all susceptibles thus continu¬ 
ously switch allegiance in the long run. In that regime. 



time 



10 10 10 10 10 
time 

FIG. 7: Typical evolution of the rescaled magnetization for 
different values of q and N with initial condition m(0) = —s 
(all B-susceptibles) on a semi-log scale: (a) Single realization 
of m{t)/s for q = 2,z = z+ = Z- = 0.2 < Zc and population 
size N = 100. At time t ~ 625, m = s and starts swinging 
back and forth the values m = ±s. Here, the MF predicts 
m*!s ~ 0.745 and ±m*/s are shown as dashed lines. Inset: 
m{t)/s vs. time for q = 3, z = 0.3 < Zc and N = 100. The 
system starts swinging between m = ±s at f « 640. Here, 
±m*/s « ±0.693 (dashed), (b) m(t)/s vs. time for the same 
parameters as in (a) but with N = 300. The system’s magne¬ 
tization switches to m = s only at f « 2 • 10®. Inset: m{t)/s 
vs. time for q = 3, z = 0.3 and N = 400. The magnetization 
switches to m = s at t ~ 3 • 10®. The difference in the switch¬ 
ing times in (a) and (b) results from the exponential scaling 
of the mean switching time on N, see text. 


the magnetization m{t) = (2n(f) — S)/N is thus charac¬ 
terized by abrupt jumps from m Ri ±s to m ~ ±s, see 
Fig.0, while the stationary ensemble-averaged magneti¬ 
zation (to(oo)) = J2m=-s '^Qm = 0 , since is even 
and each agent is as likely to be in one or the opposite 
state. A similar phenomenon has been found in the Sz- 
najd model {q = 2) with anticonformity p^ . 

This swing-state phenomenon is not captured by the 
mean field description of Sec. Ill and is here characterized 
by the mean time to switch for the first time from state 
n = 0 to n = S'. The scaling of Tq on N allows us to 
rationalize the data of Fig. ^ where the switching time is 
found to dramatically increase with the population size. 



























































FIG. 8: {Color Online) Mean switching time r® vs. z in the 
symmetric case z+ = z- = z < Zc for g = 2 (O, solid), q ~ 3 
(o, dashed) and N = 100. Symbols are results of stochastic 
simulations and lines are the predictions of dl^, see text. 


Clearly, the symmetry implies that the mean switching, 
or swinging, time Tq is identical to the mean time to 
switch from n = S to n = 0. 

Finding the mean switching time can be formulated 
as a first-passage time problem and, when N 1, Tq 
can be computed using the framework of the backward 
Fokker-Planck equation (bFPE) |^. In this context, the 
model’s bFPE infinitesimal generator is 


Gh{x) = [T+(a;) - T {x)]dx -h 


[T+{x) + T-{x)] 
2N 


dl (15) 


The mean time t^{xo) to be absorbed at a; = s (all A- 
susceptibles), starting from the initial state x = xo, with 
a reflective boundary at a; = 0 (all B-susceptibles), obeys 

0 h(xo) T^(xo) = -1, (16) 


with (d/dx)T^ (0) = 0 and t^(s) = 0 (reflective and ab¬ 
sorbing boundaries) To obtain the mean switch¬ 

ing time Tn we solve Eq. (16) with xq = 0 using standard 
methods |M, and obtain 


To" = 2N 


dy e 


-NHv) 


rv gN^iv) 
lo T+{v) + T-{v)-^ 


(17) 


where ()>(u) = -2 du As with 

other fluctuation-driven phenomena associated with 
metastable states, see e.g. |^, p7|-|^ and below, this 
result predicts that the mean switching time Tq grows 
(approximately) exponentially with the population size 
N. This explains the difference of various orders of mag¬ 
nitude in the switching time observed in Figs. IZl(a) and 

7(b). 

The predictions of (|y) are reported in Fig. g for vari¬ 
ous values oi z < Zc- These are In good agreement with 
the results of numerical simulations (averaged over 1000 
samples, each run for 10® simulation steps). When z is 
lowered well below Zc, the peaks of the SPD approach 


n = 0 and n = S. In this case, Tq Increases and switch¬ 
ing allegiance takes very long. At fixed z < we find 
that Tq increases with q. Interestingly, we also find that 
Tq can exhibit a non-monotonic dependence on z just 
below Zc when q Is kept fixed, as shown in Fig. 


B. Time-scale separation and growing fluctuations 
in the asymmetric case 


In the asymmetric case 0 < z_ < z_|_, the party A 
has more zealot supporters than party B. In this situ¬ 
ation, when z < Zc the SPD has a marked peaked near 
n = S, see Fig. ^(a). As shown in Fig. |^(b), the long time 
dynamics is characterized by a large majority of susceptl- 
bles becoming A supporters independently of the Initial 
state. The magnetization m{t) = 25z [2n{t) — S]/N 

thus fluctuates around its MF value before reach¬ 
ing m = = s -|- 2i5z when all susceptlbles are 

supporters of party A, see Fig. ||(a). The population 
composition then endlessly fluctuates, with a majority 
of susceptibles supporting party A. In this case, with 
Eq. (n), the stationary ensemble-averaged magnetiza¬ 
tion ^{oo)) = rnQ^ is positive. 

The gVM dynamics is thus characterized by various 
regimes not captured by the mean field description. For 
concreteness, we consider that the initial density of A- 
susceptibles Is xq < x*, as in Fig. ||, and distinguish four 
time scales: 

(1) After a mean time of order , the system quickly 
relaxes toward the metastable state n*_ where a random 
voter has the MF opinion ml, see Fig. ^(a). 

(ii) After a mean time rji, almost all realizations sud¬ 
denly approach the metastable state n)j_ where m{t) k, 
m^, see Figs. 11(a) and |(b). The mean transition time 
rji, as well as the average relaxation times, can be esti¬ 
mated using Kramers’ classical escape rate theory [p7| . 
The latter gives the mean transition time tk between 
the two local minima of the double-well potential U{x) 
In which an overdamped Brownian particle is moving 
subject to a zero-mean delta-correlated Gaussian white 
noise force ^{t). Here, we consider a potential U{x) such 
that dU/dx = T~{x) — T+(x), and the noise correlations 
{^{t)^{t')) =5{t-t') [T+{x*_)+T-{x*_)^. The bFPE 
generator of this Brownian particle is ( p^ with a con¬ 
stant diffusive term evaluated at x*_. Kramer’s formula 
hence gives 


T_ ~ Tk = 27r TriTra 


2N r*.* 

J X 

e 


T-(v)-T+(y) 

T-(x* ) + T+(xl) 


dy 


where Tr^ = l/y^C/"(xl), and Tra = l/\/|C/"(x*)| de¬ 
notes the mean relaxation time from state n = n* to 
n = n^. 

(iii) The system then fluctuates around before 
reaching the state n = S (all A-susceptlbles) where 
m = rnmaxi see Fig. p|(a) after a mean time r'®. In the 
realm of the bFPE, the mean time t^ for all susceptlbles 
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FIG. 9: {Color Online) (a) Typical single realization of 
mimraax = m/{s + 25z) in the asymmetric case at low 
zealotry, with q = 3 and m(0) = rrimin = — s + 2&z (all B- 
susceptibles). Here, 2 = 0.25, <5 = 0.12 and N — 1000. After a 
time t « 15 the magnetization fluctuates around ml (dashed 
line); at time t « 6T0"^ it attains mmax = s+2Sz. In principle, 
the magnetization can return to mmin (all B-susceptibles) af¬ 
ter an enormous time (-^ 10®®, not shown). Inset: Same, but 
in linear scale, (b) r® as a function of N for q = 3, z = 0.25 
and (5 = 0.12 (with m(0) = mmin). Symbols (o) are the re¬ 
sults of stochastic simulations. The solid and dashed lines are 
the predictions of (0) and respectively, showing that r’® 
grows approximately exponentially with N, see text. Inset: 
r® vs. 2 = (1005)“^ with N = 100, and Z+ — Z- — 2 kept 
fixed, for q = 2 {<>, solid) and q = 3 {o, dashed). 


to become A supporters for the first time is 


r® = 2N 


dy e' 


-NHv) 


rv ^N<t>(v) 

lo T+{v)+T-{v)- 


(18) 


When is close to the state n = S, the main contribu¬ 
tion to is given by the mean transition time that is 
independent of a:o < x*, as illustrated by Fig. ||(b). This 
is well approximated by Kramer’s formula, yielding 


r 


S 


~ 27r 


2N f®.* 

J X 

e 


T-(v}-T+(y) 
T-(x1) + T+(»!;1) 


dy 


,(19) 


showing that the mean switching time scales exponen¬ 
tially with the population size. 

(iv) The amplitude of the fluctuations around n ~ 
S, where m{t) ~ Wmax grows endlessly in time, see 


Fig. 11(b), and the system eventually returns to the state 
n = 0 (all B-susceptibles). Yet, this occurs after an enor- 

2N /®* dy 

mous amount of time, of order e (“,+)+'r (*+) ^ 

that is generally not physically observable when Af ^ 1. 

The predictions ( p^ and its approximation (|^) are 
reported in Fig. ^(b), where they are in good agreement 
with the results of stochastic simulations. These results 
confirm that r'® grows approximately exponentially with 
N when Af > 1. In Fig. |(b) , we also see that in¬ 
creases with 1 / 2 , and with q when 2 and <5 are fixed. As 
illustrated in Fig. |^(a), contrary to the case of symmet¬ 
ric zealotry, there is no “swing-state dynamics”: After a 
mean time the population persists near n ~ S where 
most susceptibles are A supporters and the magnetiza¬ 
tion is m ~ runiax, and there is virtually no switching 
back to state n ^ 0. Hence, a small bias in the zealotry, 
combined with fluctuations and nonlinearity, can greatly 
affect the voters’ opinion in the gVM. 


VI. MEAN CONSENSUS TIME IN THE 
PRESENCE OF ONE TYPE OF ZEALOTS 

When there are only A-zealots, with 2 + = C and 
2 _ = 0, an A-party consensus is always reached. Yet, the 
dynamics leading to the corresponding absorbing state 
n = S depends non-trivially on the zealotry density and 
on the initial density xq of A-susceptibles. 

Here, the fluctuation-driven dynamics is characterized 
by the mean consensus time (MCT). As illustrated in 
Fig. 1^, the MCT can change by several order of mag¬ 
nitudes when C and xq change over a small range: (i) 
Below the critical zealotry density Cc, the MCT grows 
exponentially with the population size N when xq < x*, 
see Fig. @b); (ii) Otherwise the MCT grows logarith¬ 
mically with N, see Fig. p^(inset). These phenomena are 
analyzed as follows: 

(i) When C ^ Cc and xq < x*, in line with the 
MF analysis, the density of A-susceptibles first lingers 
around the metastable state until a large fluctua¬ 
tion drives the system towards the absorbing state. This 
large-fluctuation-driven phenomenon is particularly well 
captured by the WKB theory ||2^-|^. The essence of 
this method consists of studying the quasi-stationary 
probability distribution (QSPD) 7r„ obtained by setting 
P„(t) ~ 7r„ 6"*/”'= for 0 < n < S' and Ps(t) ~ 1 - 
into the master equation (0)- The MCT is the mean de¬ 
cay time Tc of the QSPD. Since {d/dt)Ps — e“*/’’‘=/Tc ~ 
, we indeed find |^, |3^ 

A = (T+_i^s_i)-b (20) 

The computation of the MCT therefore requires finding 
the QSPD. This obeys 

+ Tn+l^n+1 ~ {Tn + = 0 , ( 21 ) 

obtained from Eq. ([|) upon neglecting an exponentially 
small term tTtlIxc- In the limit V ^ 1, the density x = 
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FIG. 10: {Color Online) (a) MCT vs. xo/s for q = 2, ^ = 0.17 
and N = 5000 (o, solid), and g = 3, ^ = 0.24 and N = 400 (o, 
dashed). Symbols are from stochastic simulations and lines 
are the predictions of (Q. (b) MCT vs. A'" for g = 2,4" = 0.17 
(o, solid), and for g = 3,^ = 0.24 (o, dashed). Initially xq = 
0.1 < a;*. Symbols (o) and (o) are from simulations (averaged 
over 10® — 10® samples), lines are the predictions of (|2^ and 
(x) are proportional to the WKB results, Eqs. ( p^ ) and (p5|). 
Inset: MCT vs. N in semi-log scale for q = 2,({ = 0.2 (o, 
solid), and for q — 3, = 0.26 and (o, dashed). 


n/N is treated as a continuous variable and Eq. ( p0| ) 
yields r"® = {tt{s)/N)\-£ r+(a:)|^^ In the continuum 


limit, Eq. (21) is solved with the WKB Ansatz 
Tr{x) - A 


( 22 ) 


where S (x) is the action, Si (x) is the amplitude, and A 

is a normalization constant [^. By substituting 
( p^ into (pl|), to leading order we find Q 


S{x)=- 'i'iy) dy, 


(23) 


where, as in Sec. IV, T'(t/) = \n[T^{y)/T~{y)]. 

Hence, when xq < x* and N 1, the leading contri¬ 
bution to the MCT is given by the accumulated action 
AS over the path joining the metastable state x = x^ 
and the unstable steady state x = x 


^N[S{x‘‘)-S(xl)] ^ gJVA5^ 


(24) 


The next-to-leading correction arising from is given 
in Refs. but for our purpose Eq. already 

provides useful information on the MCT. In fact, for q = 
2 and g = 3, the action (E^ explicitly reads 


-5(a:) = 


(1-C) 

+x In 


ln(l-C- 

A+C? 

fc(l —(■ —x) 


x) + 2( In (a: -t- C) 

(g = 2) 


C In a: -I- 2 In 
-|-(a: -b C) In 


(1-C-a;) 


(g = 3) 


(25) 


With these expressions, and with (|^) and (||) for x^ and 
a:*, the leading contribution to the WKB approximation 
of the MCT is computed explicitly, and the results re¬ 
ported in Fig[l^(b) are in excellent agreement with those 
of stochastic simulations when N ^ 1 and confirm that 
Tc grows exponentially with N. We can also check that 
AS is a decreasing function of C, which clearly implies 
that the MCT grows when ^ is decreased. 

(ii) When ^ > (jcj or for any C > 0 when the initial 
density Xq > a;*, the A-party consensus is reached much 
quicker than in the case (i), typically after a time of order 
N), see Fig. p^(inset). The backward Fokker-Planck 
formalism is again suitable to derive this result. In such 
a framework, the MCT obeys Eq. (0) supplemented by 
reflective and absorbing boundary conditions t'(0) = 0 
and t(s) = 0 Proceeding as in Sec. V, we find again 
the expression: 


,{xq) = 2N r dy [ 

Jxn 0 


V e^4>(v) dv 

0 T+{v) + T-{vy 


(26) 


As shown in the 


with (j){v) = -2/q du { r-S+rtS }- 
inset of Fig. |I^, this expression is in good agreement with 
the results of stochastic simulations and captures the 
functional dependence of the MCT whose leading contri¬ 
bution grows logarithmically with N and increases with 
q. It is also worth noting that Eq. (p^ also provides a 
meaningful approximation of the MCT in the metastable 
regime, even though when N ^ 1 its predictions are 
usually less accurate than those of the WKB method, 
see e.g. p0| . 


VII. SUMMARY AND CONCLUSION 

We have studied the dynamics of the non-linear g-voter 
model (gVM) in the presence of inflexible zealots in a fi¬ 
nite well-mixed population. In this model, voters can 
support two parties and are either “susceptibles” or “in¬ 
flexible zealots”. Susceptible voters adopt the opinion of 
a group of g > 2 neighbors if they all agree, while zealots 
are here individuals whose state never changes. The gVM 
with zealots is introduced as a simple non-trivial model 
able to capture the essence of important concepts of so¬ 
cial psychology and sociology, such as the relevance of 
conformity and independence as mechanisms for collec¬ 
tive actions and the existence of group-size thresh¬ 
old that influences the social impact of conformity |Q . 
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In spite of its simplicity and the fact that the detailed 
balance is satisfied, the dynamics of the non-linear qVM 
with zealots is rich and characterized by fluctuation- 
driven phenomena and non-trivial probability distribu¬ 
tions. The dynamics is particularly interesting at low 
level of zealotry, when the stationary distribution is bi- 
modal. In this case, we have found that when one party 
has more zealots than the other, the intensity of one peak 
greatly exceeds that of the other. The dynamics is thus 
characterized by various time scales and growing fluctu¬ 
ations around a state in which a majority of susceptibles 
support the party having more zealots. When both par¬ 
ties have the same number of zealots, below the critical 
zealotry, the long-time dynamics is characterized by the 
susceptibles endlessly swinging from a state in which they 
all support one party to the state where they all support 
the other party. We have rationalized all these features 
by computing the exact stationary probability distribu¬ 


tion and, within the backward Fokker-Planck formalism, 
the mean times for all susceptibles to switch allegiance. 
We have thus found that these mean switching times 
grow approximately exponentially with the population 
size, and they increase when the number of zealot de¬ 
creases at low zealotry. When zealots support only one 
party, we have shown that a consensus is reached in a 
mean time that grows either exponentially or logarithmi¬ 
cally with the population size, depending on the zealotry 
density and the initial condition. 

Our findings show that the properties of the nonlinear 
qVM with zealots {q > 2) are dominated by fluctuations, 
and have revealed that they are sensitive to even a small 
bias in the zealot densities. Most of the features of the 
nonlinear qVM with inflexible zealots are therefore be¬ 
yond the reach of a simple mean field analysis and gener¬ 
ally deviate from those of the classical linear voter model. 
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